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Abstract 

Atomistic simulations using a combination of classical forcefield and Density-Functional-Theory 
(DFT) show that carbon atoms remain essentially sp 2 coordinated in either bent tubes or tubes pushed by an 
atomically sharp AFM tip Subsequent Green’s-function-based transport calculations reveal that for arm- 
chair tubes there is no significant drop in conductance, while for zigzag tubes the conductance can drop by 
several orders of magnitude in AFM-pushed tubes. The effect can be attributed to simple stretching of the 
tube under tip deformation which opens up an energy gap at the Fermi surface. 
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Tremendous potential tor technological applications has thrust carbon nanotubes into one of the hottest areas 
of research activity. This has been fueled by recent experimental breakthroughs in diverse application areas 
[Articles on nanotubes in Physics World 13, Issue 6, pp. 29-53 (2000).], ranging from flat panel displays, to 
novel microelectronic devices, to hydrogen storage devices, to structural reinforcing agents, to chemical and 
electromechanical sensors. A pioneering experiment in the last application area involved a metallic nanotube 
suspended over a 600 nm long trench [T. W. Tombler et al., Nature 405, 769 (2000).]. When the middle part 
of such a suspended nanotube was pushed with the tip of an Atomic Force Microscope (AFM), the conduc- 
tivity was found to decrease by nearly two orders of magnitude for a deformation angle of 15°. The effect 
was found to be completely reversible, i.e., through repeated cycles of AFM-de formation and tip removal, 
the electrical conductance displayed a cyclical variation with constant amplitude. 

The drop in conductance in the AFM-deformed tube was much higher than the computationally 
predicted values for tube; bent under mechanical duress. Such calculations, using both tight-binding [M. 
Nardelli and J. Bemholc. Phys. Rev. B 60, R 16338 (1999).] and semi-empirical Extended-Htickel type 
approaches [A. Rochefort P. Avouris, F. Lesage, and D. Salahub, Phys. Rev. B 60, 13824 (1999).] concluded 
that even under large bending angles the reduction in electrical conductance was less than an order of magni- 
tude. For AFM-deformed nanotubes, in contrast, O(N) tight-binding calculations [L. Liu et al., Phys. Rev. 
Lett. 84, 4950 (2000).] show that beyond a critical deformation several C-atoms close to the AFM tip 
become sp 3 -coordinated. The sp 3 coordination ties up delocalized n-electrons into localized a-states. This 
would naturally explain the large drop in electrical conductivity, as verified by explicit transport calcula- 
tions. 

Under either bending or pushing by an atomically sharp AFM tip, bond reconstruction, if any, is 
likely to occur only in the highly deformed, non-straight part of the tube in the middle. This prompted us to 
use a DFT-based quantum mechanical description of the middle part of the tube (~ 100 atoms), while the 
long and essentially straight part away from the middle was described accurately using the Universal Force- 
field (UFF) [A. K. Rappe et al., J. Am. Chem. Soc. 1 14, 10024 (1992)., N. Yao and V. Lordi, J. Appl. Phys. 
84, 1939 (1998).]. Structures and energetics obtained this way for bent tubes were in good agreement with 
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previous work using an interatomic potential [S. Iijima. C. Brabec. A. Main, and J. Bemholc, J. Chem. Phys. 
104. 2089 { 1996).]. For the AFM tip-deformed tubes, on the other hand, the situation depended on how the 
tip was represented. Thus, if the presence of the tip was simulated by constraining a single C-atom on the 
bottom side of the middle part of the tube, simulations on a (5, 5) armchair led to the development of sp 3 
coordination between the constrained atom and an atom on the top side at a critical deformation angle of - 
7°, which destabilized into a complex broken-bond defect at higher deformation angles [A. Maiti, Chem. 
Phys. Lett. 331, 21 (2000).]. On the other hand, a more realistic representation of the AFM-tip by means of 
an atomically sharp 15-atorn Li-needle yielded an sp 2 -coordinated all-hexagonal tube for deformation angles 
as high as 25* for the same (5, 5) tube [A. Maiti, Phys. Stat. Sol. B 226, 87 (2001).]. A large drop in conduc- 
tance is expected under sp* coordination and broken bond defects [L. Liu et al., Phys. Rev. Lett. 84, 4950 
(2000).]. However, given the uncertainty of sp^ coordination, can one still expect a significant conductance 
drop in a tip-deformed, yet sp~ -coordinated tube ? 

In this Letter, we address the above question by extending the combined DFT-UFF calculations 
to a (12, 0) metallic zigzag tube. For comparison, we have also considered a (6, 6) armchair tube, which is 
slightly smaller in diameter, but has the same number of atoms ( 1 2) along the tube circumference as the ( 1 2, 
0) zigzag. The main result of structural relaxation using DFT-UFF is that, sp 3 coordination does not happen 
under either bending or tip- deformation (using an atomically sharp 15-atom Li-tip as in ref. [A. Maiti, Phys. 
Stat. Sol. B 226, 87 (2001 ).]) up to very large angles. Following structural relaxation at each bending and 
deformation angle, we compute the electronic density of states (DOS), transmission and conductance using 
the recursive Green’s function method [A. Svizhenko, M. P. Anantram, T. R. Govindan and B. Biegel, Pre- 
print, Submitted to JAP.]. For the armchair tube, the resulting conductance is lowered only by a factor of 
1.0 1 for the tube bent by 40° and a factor of 1.05 for a tube tip-deformed by 25°. Under the same deforma- 
tions, drop in conductance for the zigzag tube, is much higher, being a factor of 1.9 under bending, and a 
remarkable 1.7-1 0 4 under t;p-deformation. 

The simulations were carried out on tubes of 2400 atoms, both for the (6, 6) and (12, 0) tubes. 
Initially the straight tube v/as relaxed with the UFF. For bending simulations, two halves of the tube were 
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then rotated by equal and opposite angles about an axis perpendicular to the tube and passing through the 
center of mass of the initial straight tube. At each end of the tube, a contact region defined by a unit cell [For 
computation of self-energies of semi-infinite carbon nanotube contacts in Eqn. ( l ), it was convenient to par- 
tition the whole tube in adjacent repeating segments, or “unit cells”. For an armchair and a zigzag tube 
respectively, the unit cell consists of 2 and 4 rings of atoms around the circumference. This implies 24 atoms 
for the (6, 6) tube and 48 atoms for the (12, 0) tube per unit cell.] plus one atomic ring (a total of 36 and 60 
atoms for the armchair and the zigzag tube respectively) was then fixed and the whole tube relaxed with the 
UFF. To simulate AFM-tip-deformation, the 15-atom Li-needle was initially aimed at the center of a hexa- 
gon on the bottom-side of the middle part of tube. The Li-needle was then displaced by an amount 5 toward 
the tube along the needle-axis, resulting in a deformation angle 9 = tan' 1 (28/L), L being the unstretched 
length of the tube. The whole tube was then relaxed by UFF keeping the needle atoms and the end contact 
regions of the tube fixed. Fixing the relative positions of contact region atoms at the same value as in an 
unstretched tube guarantees that contacts may be approximated by ideal undeformed semi-infinite carbon 
nanotube leads and that all possible contact modes are coupled to the deformed part of the tube. 

Following the UFF relaxation, a cluster of 132 C-atoms for the (6, 6) and a cluster of 144 C- 
atoms for the ( 12, 0) were out out from the middle of the tubes. These clusters, referred to below as the QM 
clusters (plus the 15 Li-tip atoms in tip-deformation simulations) were further relaxed with Accelrys’ DFT- 
code DMol 3 [Delley, B. J. Chem. Phys. 1990, 92, 508; J. Phys. Chem. 1996, 100, 6107; http://www.acce!- 
rys.com/mstudio/dmol3.httnl.], with the end atoms of the cluster plus the Li-tip atoms fixed at their respec- 
tive classical positions. The electronic wave functions in DMol 3 were expanded in a double-numeric 
polarized (DNP) basis set with a real-space cutoff of 4.0 A. The Hamiltonian was approximated with the 
Harris functional [Z. Lin and J. Harris, J. Phys.: Condens. Matter 4, 1055 (1992).] using a local exchange- 
correlation potential [S. H. Vosko, L. Wilk, and M. Nusair, Can. J. Phys. 58, 1200 (1980).], and the 
“medium” grid was chosen for numerical integration. 

Fig. 1 displays the tip-deformed QM-cluster for the (6, 6) and the (12, 0) tubes at the highest 
deformation angle of 25° considered in these simulations. Even under such large deformations, there is no 
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indication of sp J bonding, similar to what was observed for the (5. 5) tube in ref. [A. Maiti. Phys. Stat. Sol. 
B 226. 87 (2001).]. Although not explicitly shown here, results for bending also yield sp : -coordinated all- 
hexagonal tubes. The absence of sp 3 coordination is inferred based on an analysis of nearest neighbor dis- 
tances of the atoms with the highest displacements, i.e., the ones on the top of the kink in a bent tube, and the 
ones closest to the Li- tip in a tip-deformed tube. Although for each of these atoms the three nearest neighbor 
C-C bonds are stretched to between 1.45-1.75 A, the distance of the fourth neighbor, required to induce sp 3 
coordination is greater than 2.2 A for all tubes in our simulations. The main difference between a tip- 
deformed tube versus a bent tube is that there is an overall stretching in the former [For a tube with a very 
large length-to-diameter ratio, the length L stretches to - L sec q, q being the tip-deformation angle. How- 
ever, for moderately long tubes used in our simulations, the average tensile strain in the straight part of the 
tube is slightly lower than ( sec q - 1 ).] whereas in the latter case there is no net stretching, and the extra com- 
pressive strain on the bottom side is relieved through the formation of a kink beyond a critical bending angle. 

Following atomic relaxation of the structures, we performed conductance calculations in order 
to make further predictions on the electromechanical behavior of nanotubes. A coherent conductance was 
studied within a nearest-neighbor sp 3 -tight-binding Hamiltonian in a non-orthogona! basis. The parameter- 
ization scheme explicitly accounts for effects of strain in the system through a bond-length-dependence of 
the Hamiltonian and the overlap matrices // J; and 5^, as in Ref [D. A. Papaconstantopoulos, M. J. Mehl, S. C. 
Erwin and M. R. Pederson Tight-Binding Approach to Computational Materials Science, P.E.A. Turchi, A. 
Gonis, and L. Colombo, eds., MRS Proceedings 491, (Materials Research Society, Warrendale, PA, 1998)]. 
We have also checked to confirm that other tight binding parameterizations give qualitatively the same 
results [J.-C. Charlier, Ph. Lambin and T. W. Ebbesen, Phys. Rev. B 54, R8377 (1996)., W. A. Harrison, 
Electronic Structure and the Properties of Solids (Dover, New York, 1989)]. First, the retarded Green's func- 
tion C R of the whole nanotube was determined by solving the following equation: 

(E-S, l -H, l -I. Ljl -I. Rii )G Rjk =8,‘ t ([) 

where are the retarded self-energies of the left and the right semi-infinite contacts. The transmission 

and the electronic density of states (DOS) at each energy were then found [In the Eqs.(l-3), summation is 
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performed over the repeating roman indices. The lower and upper indices denote covariant and contravariant 
components of a tensor., D Lohez and M. Lanoo, Phys. Rev. B 27, 5007 (1983 ).| from the equations: 


nE) = G l<JI Y Ljk G AM Y KM , (2) 

N a (E) = — 

* , ( 3 ) 


where T^# =i(L R ^ ^ ) are the couplings to the left and right leads. Finally, the total conductance of 

the tube was computed using Landauer-Biittiker formula: 


2e 2 7 3/ 

G ~ [WW-fedE 


( 4 ) 


where fjE) is the Fermi-Dsrac function. 

Fig. 2 displays the computed conductance (at T = 300K) for the (6, 6) and the (12, 0) tubes as a 
function of bending and tip-deformation angles. The conductance remains essentially constant for the (6, 6) 
tube in either bending or tip-deformation simulations. However, for the (12, 0) tube the conductance drops 
by a factor of 1.9 under bending at 9=40°, and much more significantly under tip deformation: by - 0.3 at 
15°, two orders of magnitude at 20°, and 4 orders of magnitude at 0=25’ [For the (12, 0) tube, we also stud- 
ied the dependence of conductance-drop as a function of tube-length. Thus at a tip-deformation angle of 15 
o, tubes of 2400, 3600 and 4800 atoms have conductance-drops of 0.29, 0.13 and 0.10 respectively, which 
extrapolates to - 0.08 for very long tubes.]. 

To analyze which part of the zigzag tube is responsible for the conductance drop, we computed 
the DOS in the vicinity of the Fermi surface. Fig. 3 displays the DOS averaged over 2 unit cells [For compu- 
tation of self-energies of semi-infinite carbon nanotube contacts in Eqn. (1), it was convenient to partition 
the whole tube in adjacent repeating segments, or “unit ceils 1 ’. For an armchair and a zigzag tube respec- 
tively, the unit cell consists of 2 and 4 rings of atoms around the circumference. This implies 24 atoms for 
the (6, 6) tube and 48 atoms for the (12, 0) tube per unit cell.] (96 atoms) in three different regions of the 
AFM-deformed ( 1 2, 0) tube for 0=25°: ( 1 ) undeformed contact, (2) highly deformed tip region, and (3) the 

uniformly stretched straight regions on either side of the tip-deformed region. DOS in both tip region and 
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the straight part show a bandgap opening, which proves that the conductance drop occurs everywhere tn the 
tube, rather than in the tip-deformed region alone. This has important implications tor the application of nan- 
otubes as electromechanical sensors: given a metallic zigzag nanotube, one could induce a significant con - 
ductance drop simply upon uniform stretching. To check this, we computed the conductance of a uniformly 
stretched (12, 0) tube as a function of strain, shown in Fig. 4, and compared it to that of AFM-deformed tube 
from Fig. 2. Both cases show quantitatively the same drastic decay of conductance. The inset in Fig. 4 also 
shows a band gap opening in transmission in the two cases [A. Heyd, A. Charlier and E. McRae, Phys. Rev. 
B 55, 6820 (1997).], compared to that of a non-deformed tube. 


In order to explain the differences in conductance of the (6, 6) and the ( 1 2, 0) tubes as a function 
of strain, we have analyzed the bandstructures of a metallic zigzag and an armchair nanotube. Starting from 
the band structure of the 2D graphene [J. W. Mintmire and C. White, Carbon, V.33, No. 7. pp. 893-902 
(Elsevier Science Ltd., 1995); M. S. Dresselhaus, G. Dresselhaus, and P. C. EkJund, Science of Fullerenes 
and Carbon Nanotubes (Academic Press, New York, 1996).] under deformation [For electromechanical 
properties of tubes of any chirality and definition of tl, t2, t3, see, L. Yang, M. P. Anantram, J. Han and J. R 
Lu, Phys. Rev. B 60, 1387^ (1999).], one can derive the following dispersion relations for the crossing sub- 
bands within the n-electron approximation (a' below is the strained periodic repeat length along the nano- 
tube axis): 


E(k) = ±t 2 {\ 


■> n 3kci . \f-i 

+cr -2acos( )} 


(5a) 


for a metallic zigzag tube, and 

E(k) = ±r. 


|1 — 2ct cos(— ) 
2 


(5b) 


for an armchair tube. Here, a=t/t 2 is the ratio of nonequivalent hopping parameters. For a zigzag tube, a>/. 
Consequently, there is no value of ka ' for which E(k) = 0 (for dispersion (5a)), and a band gap opens up. 
However, for an armchair tube, a< 1 and one can always find a value of ka \ such that E(k)=0 (for dispersion 
(5b)), as long as a>//2. The magnitude of the strain-induced bandgap decreases monotonically with 


increase in chiral angle [For electromechanical properties of tubes of any chirality and definition of tl, t2, t3, 
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see. L. Yang, M. P. Anartram. I. Han and J. P. Lu. Phys. Rev. B 60. 13874 (1999).], being maximum for a 
chiral angle of 0° (zigzag) and gradually reducing to zero at a chiral angle of 30° (armchair). An experiment 
as in ref. [T. W. Tombler «t aL Nature 405, 769 (2000).] is, thus, expected to show a decrease in conductance 
as the nanotube is deformed with an AFM tip, for all nanotubes except the armchair tube. The decrease in 
conductance should be large for the metallic zigzag nanotubes, and small for nanotubes with a chiral angle 
closer to armchair. 

In summary, we find that both under bending and under deformation with an atomically sharp 
AFM-tip, carbon nanotubes essentially remain all-hexagonal and sp~-coordinated. In the absence of sp 3 
coordination, armchair tubes remain significantly conducting even at large deformations. However, metallic 
zigzag tubes display a dramatic drop in conductance, particularly under tip-deformation. A density of states 
analysis indicates that the conductance drop is distributed over the whole tube, rather than focused in the tip 
region. This suggests the possibility of designing nanoelectromechanical sensors in which nanotubes are 
subjected to a uniform tensile strain. 
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Figure captions : 


Fig I. DMol J -reIaxed Li-tip -deformed QM clusters for: (a) the {6. 6) armchair ( 132 C-atoms); and (b) the 
( 1 2, 0) zigzag ( 144 C-atoms). in side views. The deformation angle is 25® for both tubes. Figs, (c) and (d) are 
respective views along the cube length, with the Li-tip hidden for clarity. 

Fig2. Computed conductance (at T = 300K) for the (6, 6) and (12, 0) nanotubes as a function of: (top) bend- 
ing; and (bottom) tip-deformation. Under 40° bending conductance of the zigzag tube drops by a factor of 
1.9, while for the armchair tube it drops by only 1 .01. Under 25° tip-deformation, conductance of the zigzag 
tube drops by 4 orders of magnitude, while for the armchair tube it drops only by a factor of 1 .05. 

Fig. 3. Density of states, averaged over two unit cells (96 atoms) in three different regions of AFM-deformed 
(12. 0) tube at 0 = 25°. Thr Fermi surface is at E=0. Nearly equal band gaps open up in the straight part of 
the tube and in the tip region. Lower DOS in the tip region is due to a larger local straining of C-C bonds 
close to the AFM-tip. 

Fig. 4. Conductance of the uniformly stretched (12, 0) tube compared to that of the tip-deformation case in 
Fig. 2. Actual angles of tip-deformation are indicated. The % strain for the AFM-deformed tube is com- 
puted from the average C-C bond-stretch in the middle of the straight portion of the tube [For a tube with a 
very large length-to-diame^er ratio, the length L stretches to - L sec q, q being the tip-deformation angle. 
However, for moderately long tubes used in our simulations, the average tensile strain in the straight part of 
the tube is slightly lower than (sec q - I).]. The inset shows transmission in the vicinity of the Fermi surface 
(E = 0) for a uniform strain of 10% and a tip-deformation angle of 25°, as compared to the undeformed tube. 
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